> 

(N 

a^ 
o 
o 



^ 



Oh! 



X 



EIGENVALUES OF THE QCD DIRAC OPERATOR AT FINITE 
TEMPERATURE AND DENSITY 



E. BITTNER'', M.-P. LOMBARDO'', H. MARKUM^', R. PULLIRSCH'^ 



O 

A—^ I ^Institut fiir Kernphysik, Technische Universitat Wien, 

fvj ' A-IO4O Vienna, Austria 

Istituto Nazionale di Fisica Nucleare, 



Sezione di Padova, e Gruppo Collegato di Trento, Italy 



Oh' 
(D 

^-* ' We investigate the eigenvalue spectrum of the staggered Dirac matrix in two-color 

|/~\ ' QCD at nonzero temperature and at baryon density when the eigenvalues become 

complex. The quasi-zero modes and their role for chiral symmetry breaking and the 
deconfinement transition are examined. The bulk of the spectrum and its relation 
to quantum chaos is considered. Comparison with predictions from random matrix 
theory is presented. 



1 Chiral Condensate 

The properties of the eigenvalues of the Dirac operator are of great in- 
terest for important features of QCD. The Banks-Casher formulaEl relates 



f^ ■ the Dirac eigenvalue density p(A) at A = to the chiral condensate, S 



KV'V')! — hnig^o limy^oo 7rp(e)/l^. The microscopic spectral density, Ps{z) = 

O ; limy_>oo P (z/VT,) /V^ should be given by the appropriate result of randoin 

matrix theory (RMT),I3 which also generates the Leutwyler-Smilga sum rules. 

A formulation of the QCD Dirac operator ^at chemical potential /.t ^ on 

the lattice in the staggered scheme is given byQ 



Mx^y{U,fi) = — ^ [U^{x)t]i, {x)Sy ^x+i, -h.c] 

+ ^ [Ui{x)vi{x)e'^Sy^^^i - Ul{y)f^i{y)e-'^6y^,_i\ , (1) 

with the link variables U and the staggered phases 77. We report on computa- 
tions with gauge group SU(2) on a 6^ lattice a.t (3 — 4/^^ = 1.3 and with Nf — 2 
flavors of staggered fermions of mass m — 0.07. For this system the fermion 
determinant is real and lattice simulations of the full theory with chemical po- 
tential become feasible exhibiting a phase transition at fie ~ m.^/2 k, 0.3 where 
the chiral condensate (nearly) vanishes and a diquark condensate develops. EI 
In the left plot of Fig. 1 we compare the densities of the small eigenvalues 
at /i = to 0.4 on our 6'* lattice, averaged over at least 160 configurations. 
Since the eigenvalues move into the complex plane for /i > 0, a band of width 
e = 0.015 parallel to the imaginary axis is considered to construct p{y), i.e. 
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Figure 1: Left plot: Density p(y) of small eigenvalues for two-color QCD on a 6* lattice 
from fi = to 0.4. The loss of quasi-zero modes is accompanied by a vanishing of the chiral 
condensate. Right plot: Chiral condensate extracted by three different methods (see text). 



p{y) = J_ dx p{x,y), where p{x,y) is the density of the complex eigenvalues 
X + iy. 

The density p{y) is used to estimate a value for the chiral condensate by 
naively applying the Banks-Casher relation which originally was derived for 
real eigenvalues appearing in pairs of opposite sign. We further employed the 
standard definition of the Green's function^ by inverting the fermionic matrix 
with a noisy source and by computing its eigenvalues exactly, respectively, to 
get the condensate. Thus the chiral condensate for two-color QCD with finite 
chemical potential was extracted by three methods. The preliminary results 
for < ^tjj > and its variance are shown in the righthand plot of Fig. 1. 



2 Quantum Chaos 

The fluctuation properties of the eigenvalues in the bulk of the spectrum have 
also attracted attention. It was shown in Ref. g for Hermitian Dirac opera- 
tors that on the scale of the mean level spacing they are described by RMT. 
For example, the nearest-neighbor spacing distribution P{s), i.e. the distribu- 
tion of spacings s between adjacent eigenvalues on the unfolded scale, agrees 
with the Wigner surmise of RMT. According to the Bohigas-Giannoni-Schmit 
conjecture, lI quantum systems whose classical counterparts are chaotic have a 
nearest-neighbor spacing distribution given by RMT whereas systems whose 
classical counterparts are integrable obey a Poisson distribution, -Pp(s) = e~*. 
Therefore, the specific form of P{s) is often taken as a criterion for the presence 
or absence of "quantum chaos" . 

For ^ > 0, the Dirac operator loses its Hermiticity properties so that its 
eigenvalues become complex. The aim of the present analysis is to investigate 
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Figure 2: Nearest-neighbor spacing distribution for two-color QCD with varying /i. The 
analytic curves are the Wigner distribution, Pw = 262144/(7297r^)s'' exp(-64/(97rls^) (left), 
the Ginibre distribution of Eq. (Eh (middle) and the Poisson distribution of Eq. (pi) (right). 



whether non-Hermitian RMT is able to describe the fluctuation properties 
of the complex eigenvalues of the QCD Dirac operator. We apply a two- 
dimensional unfolding procedure^ to separate the average eigenvalue density 
from the fluctuations and construct the nearest-neighbor spacing distribution, 
P{s), of adjacent eigenvalues in the complex plane. Adjacent eigenvalues are 
defined to be the pairs for which the Euclidean distance in the complex plane is 
smallest. rThe data are then compared to analytical predictions of the Ginibre 
ensembleu of non-Hermitian RMT, which describes the situation where the 
real and imaginary parts of the strongly correlated eigenvalues have approx- 
imately the same average magnitude. In the Ginibre ensemble, the average 
spectral density is already constant inside a circle apd zero outside. In this 
case, unfolding is not necessary, and P(s) is given byllj 



Pg{s) — cp{cs) , p{s) ~ 2s lim 
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where e„(x) — X]m=o^™/"^' ^-'^'^ '^ ~ /o dssp{s) — 1.1429.... For uncoixe- 
lated eigenvalues in the complex plane, the Poisson distribution becomesE£l 



Pp(s) 



-rrs^/4 



(3) 



This should not be confused with the Wigner distribution for a Hermitian 
operator. Q 

Our results for P{s) are presented in Fig. 2. As a function of /i, we expect 
to find a transition from Wigner to Ginibre behavior in P{s). This was clear Ik 
seen in color-SU(3) with Nf = S flavors and quenched chemical potential, H 



where differences between both, curves are more pronounced. For the symplec- 
tic ensemble of color-SU(2) with staggered fermions, the Wigner and Ginibre 
distributions are very close to each other and thus harder to distinguish. They 
are reproduced by our preliminary data for /i == and fi ~ 0.4, respectively. 
Even in the deconfined phase, where the effect of the chemical potential might 
order the system, the gauge fields retain a considerable degree of randomness, 
which apparently gives rise to quantum chaos. 

For ^ > 1.0, the lattice results for P{s) deviate substantially from the Gini- 
bre distribution and could be interpreted as Poisson behavior, corresponding 
to uncorrelated eigenvalues. A plausible explanation of the transition to Pois- 
son behavior is provided by the following two (related) observations. First, for 
large /x the terms containing e^ in Eq. pj) dominate the Dirac matrix giving 
rise to uncorrelated eigenvalues. Second, for large n the fermion density on 
the finite lattice reaches saturation due to the limited box size and the Pauli 
exclusion principle. 
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